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Abstract 

We present a first calculation of the dilepton yield and elliptic flow done with 3+1D viscous 
hydrodynamical simulations of relativistic heavy ion collisions at the top RHIC energy. A com- 
parison with recent experimental data from the STAR collaboration is made. 

Si 

Q_! 1. Introduction 

^-l The lepton pairs resulting from nuclear collisions at energies typical of the Relativistic Heavy 

(— i | Ion Collider (RHIC) and of the Large Hadron Collider (LHC) will come from various sources 

which include high-temperature QCD processes and reactions involving thermal hadrons in the 
confined sector. The decay of open-charm mesons constitutes an irreducible background which 
must be included in theoretical analyses, in the absence of direct vertex detection. The analyses 
now available from RHIC and soon to be done at the LHC will provide a wealth of observables to 
qualify the enduring success of relativistic hydrodynamics, and to verify whether this successful 
phenomenology extends to electromagnetic observables. We present the results of a study of 
dilepton production from a hot environment with finite shear viscosity. 

2. Thermal dilepton production rates and their viscous correction 

In this work, the quark-antiquark annihilation at leading order (the Born approximation) 
is used as the dilepton production rate at high temperatures. To include shear viscosity, the 
thermal distribution functions n are augmented to include a correction 5n{p) — C- s jpn{p){\ + 
n(p))p a pP^f L d]. The rates now depend on 7r /JV , the viscous correction to the stress-energy tensor, 
and therefore on details of the dynamical viscous hydrodynamical simulations. We use music, the 
three-dimensional hydrodynamical simulation of Ref. yjj. In the phase where composite hadrons 
exist and interact, the dilepton thermal production rate is related to the in-medium behaviour of 
the current-current correlator J3l, the Wightman function. Using the Kubo-Martin-Schwinger 
(KMS) relation, one can write the lepton pair production rate, R, in terms of the imaginary part 
of the retarded photon self-energy, 
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where the lepton rest masses have been set to zero, a is the electromagnetic fine structure con- 
stant, and q~ = M' , M being the dilepton invariant mass. An equation involving the imaginary 
part of the retarded in-medium vector propagator, Im D^, is obtained from Eq. ([TJ using the Vec- 
tor Dominance Model (VDM). Here, we will not extend the definition of the Wightman function 
nor that of the KMS relation in the presence of viscosity. Instead we explore the changes viscos- 
ity entails by modifying the thermal self-energy in the propagator and the Bose distribution in 
Eq. flj, both via 6n. The total self-energy is IT"' = II™ C + U} + 6T§. The calculations of n™° 
are presented in Ref. |3.|5l0]; there, effective Lagrangians describe all interactions contributing 
to II™ C . For the in-medium self-energy Tl v , we are using here the approach in Refs. , where 
a forward scattering amplitude fvais) is used to obtain the thermal self-energy of the vector me- 
son V. The self-energies of Ref. |8[] will be considered in upcoming work. Finally, dYly is the 
correction to the self-energies that stems from the inclusion of viscous effects. 

3. Dilepton production from charm decays 

In proton-proton and Au+Au collisions with center-of-mass energies of 200 GeV, dileptons 
produced from the semi-leptonic decay of pairs of charm quarks dominate the yield in the "in- 
termediate mass range" (between the and the J /(//). We use pythia8 to generate events with 
heavy quarks. Note that pythia includes the important processes of flavour excitation as well as 
the radiative splittings important for large pj. We also use eks98 to determine the initial par- 
ton distribution functions in the nuclei. Then, using the same hydrodynamical description as was 
used to determine the thermal dilepton production, the heavy quarks are evolved using relativistic 
Langevin dynamics and the heavy quark spatial diffusion coefficient D c = 3/{2jtT). The heavy 
quarks then hadronize. The procedure is given in more details in Ref. J^t]. 



4. Results 
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Figure 1: (a) Dilepton yield from the the hadronic and QGP phases as a function of invariant mass, (b) Dilepton yield 
from the QGP extending to high invariant masses: ideal (solid curve) and viscous (dashed curve) are presented. All 
results shown here are for a minimum bias centrality class. 



We show results with viscous and inviscid hydrodynamics, with to = 0.4 fm/c, and consider 
Au + Au collisions at the top RHIC energy: V^nn = 200 GeV. In the viscous simulations, the 
generation of entropy will impose a lower initial temperature in order to reproduce the final state 
hadronic observables. The ratio of shear viscosity to entropy density is set to rj/s = l/An. 
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Figure Q] (a) presents the yield of dielectrons from the HG and QGP phases, with viscous 
effects. For the viscous HG calculation, we verified that the viscous corrections have little to no 
effect on the dilepton yield by correcting the numerator and denominator of Eq. (0 together, and 
individually. This result is in line with that obtained for real photons ifioll . and it is consistent 
with a small n MV in the HG phase. In the QGP phase, where if x is more significant, the shape 
should be changed in the region where the viscous correction to the thermal rate is large (at large 
invariant masses). This can be seen for M ~ 2 - 3 GeV in Fig. [T](b). However, shear viscosity 
effects on dilepton yields remain modest. 

The elliptic flow of dileptons is more sensitive to viscous effects than their spectra, as was 
also found for real photons ifioll . Figure [2] shows that the peaks around the p, a»,and <p masses 
first seen in a calculation using vacuum spectral distributions 1.1 ill are present in viiM), but are 
significantly reduced owing to in medium effects. The shear viscosity effects on the dynamics 
slightly broadens the M distribution (Fig. 2 (b)), as the temperature drop there is slower ifioll . 
Note finally here that our \'2(M) calculation is approximate, in the absence of a formal extension 
of the KMS relation to out-of-equilibrium environments. 




5. Comparison with STAR data 

The plots in Figure [3] show the dielectron yield and \>2 as functions of invariant mass, now 
including the contribution from charm decays. Heavy quark energy loss has a significant effect 
on the dilepton yields in the intermediate mass range: the invariant mass spectrum is reduced 
by almost an order of magnitude. In addition, the flow of heavy quarks leads to an additional 
azimuthal anisotropy of the lepton pairs. Including the effect of energy loss leads to a better 



agreement with the current data on dN/dMdy for intermediate mass dileptons from STAR 11211 . 
It is fair to write, however, that both scenarios are still consistent with the data, given the current 
experimental error bars. Finally, not included in this first analysis are the possible effect of higher 



mass thermal hadrons on the yield and on the elliptic flow of intermediate mass dileptons II 1311 



6. Conclusion 



We have presented a calculation of lepton pair production made using a 3+ ID relativistic 
hydrodynamics approach, with finite shear viscosity. We have shown the effect of charmed 
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Figure 3: (a) Comparison of the dilepton yield as a function of M with experimental data from STAR fT3l . (b) Dilepton 
i>2 versus invariant mass including charmed hadrons. All data and calculations shown here are for a 0-10% centrality 
class. 



quark energy loss on dilepton spectra and elliptic flow. We compared with recent STAR data. 
Future efforts will explore alternate vector spectral densities, will exploit recent advances in our 
knowledge of initial states, and will consider LHC conditions. 
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